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Abstract 

The particle approach to one-dimensional potential scattering is applied to non relativistic tun- 
nelling between two, three and four identical barriers. We demonstrate as expected that the 
infinite sum of particle contributions yield the plane wave results. In particular, the existence of 
resonance/transparency for twin tunnelling in the wave limit is immediately obvious. The known 
resonances for three and four barriers are also derived. The transition from the wave limit to the 
particle limit is exhibit numerically. 

I. INTRODUCTION. 

The particle approach to piece-wise potential scattering is characterized by considering the reflection 
and transmission amplitudes, successively, at each potential discontinuity, this method, applied also 
in optics [1], can be used as an alternative to plane vifave continuity equations for the derivation 
of the transmission and reflection amplitudes and find application in approximate solutions to one- 
dimensional potential problems [2,3]. Several interesting papers have analyzed tunneling vi^itli this 
method [4-6]. However, in certain situations, this approach is the natural physical choice in the 
"particle limit" in which the incoming wave packets are small compared to the potentials extension. 
However, it must be emphasized that the approach itself does not require the use of wave packets. It 
predicts multiple reflections and even for a single incoming wave it generally results in infinite reflected 
and transmitted waves. An example of this approach is the diffusion above a single potential barrier, 
from which one significant consequence, otherwise mysterious, is the step limit for large barrier [7]. 
The first barrier reflection coefficient reproduces the step result. 

The alternative and standard approach to potential scattering is with a single wave analysis [8]. 
We refer to this as the "wave limit" . It is characterized by continuity equations often best described 
by the use of matrices and yields a single reflected and transmitted amplitude. Due to the fact that 
the algebraic sum of the infinite contributions in the particle approach yields the wave result, the two 
methods are mathematically equivalent. However, they are not equivalent in practice because a sum 
of calculated terms, if finite, is an unambiguous process, but the decomposition of an expression as 
an infinite sum, on the other hand, is highly ambiguous, even if only one decomposition can lead to 
probability conservation and to the correct exit times of the separate wave packets. In the wave limit, 
probabilities arc calculated by first summing the particle terms and then squaring the modulus. In 
the particle limit, the probabilities arc calculated by first squaring the modulus of each term and then 
summing. 

Resonance phenomena, within the Schrodinger equation, are well known for a single barrier when 
considering plane waves with energy above the potential value, i.e. E > Vq. It results in unit 
transmission probabilities. For a given plane wave energy there are unlimited such resonances as the 
barrier length increases. More realistically, for an incoming wave packet it can be shown that resonance 
occurs only if the barrier length is much smaller than the wave packet dimensions [9] . Resonance is 
a property of the wave nature of particles while, in the other limit in which the wave packet is small 
compared to the barrier length, the transmission amplitude breaks up into infinite transmitted (and 
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consequently reflected) wave packets, the so called "particle limit". There is no particle limit for 
tunnelling, E < Vq, through a single barrier [10-12]. However, for more than one barrier there exists 
another coherence effect, tunnelling resonances [13-18]. For particular values of the momentum and 
of the inter-barrier distance, identical barriers become transparent, i.e. the tunnelling transmission 
probability equals unity. 

In the next section, we recall the single barrier tunnelling results and define the amplitudes relevant 
to the particle approach. We then consider the case of twin barriers. We illustrate our procedure 
for calculating the individual particle contributions and then sum them to obtain the wave limit 
result. This sum can also be viewed as a compact expression for the particle series. The resonance 
condition will be obvious in this expression. We shall also prove probability conservation in both the 
wave and particle limit. Conservation of probability must also be valid for all intermediate (partial 
wave packet overlap) cases, however no simple analytic calculation of these probabilities is known. 
In section III, we describe some numerical calculations in which an incoming Gaussian wave packet 
is used. The existence, albeit degraded, of resonances is indicted by oscillations as a function of 
inter-barrier distance d. For large d, the particle picture sets in and the oscillatory behavior is damped 
out. In section IV, we consider the three identical barrier case and derive the, non obvious, resonance 
conditions. To proceed to higher numbers of barriers it is more practical to change the method of 
calculation. In section V, we derive a matrix expression for the wave limit and then describe how this 
can be rewritten in a form from which the individual particle series can be extracted. This is applied, 
as an example, to the case of four identical barriers and the first transmitted particle terms confirm 
the direct particle procedure. We draw our conclusions in section VI. 



II. TWIN BARRIER TUNNELLING 



To perform our calculation, we need the refiection and transmission amplitudes for a single barrier. 
The single barrier reflection amplitude depends not only upon the barrier height and length but also 
on the barrier position and consequently upon whether the barrier is encountered coming from the 
left or from the right. Nevertheless, the explicit rules are quite simple. 

Let R be the reflection amplitude for a plane wave of unit amplitude with energy E < Vo impinging 
on a barrier of length L from the left at x ~ 0, see Fig. 1-a. Let T be the transmission amplitude for 
the same wave. Notice that the refiection and transmission probabilities, also known as the reflection 
and transmission coefficients, are consequently \R\ and \T\ respectively. Now, a standard plane wave 
calculation gives the results 
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R = —i COS0 tanh(pL) exp(i0) = — i |i?| exp(i0) , (1) 

2kp 

^ = ^?^cxp[*(</)-fcL)] = |r|exp[*(</)-fcL)] , (2) 

with 



k = V2mE/n , p = ^/2m{Vo - E) /h and tan = (fc" - p^) tanh(/9L)/ 2fcp , 
where, by convention, 

-7r/2 < (/)< 7r/2 . 

Conservation of probability results from the fact that \R\ + \T\ =1. Now, shift the front of the 
barrier from the origin to the position x = a. The reflection amplitude acquires an additional phase, 
becoming 

R exp(2 ik a) . (3) 
If the wave impinges the barrier from the right (momentum — fc), we obtain instead 

R cxp{-2 ikb) , (4) 

where b = a + L. In both cases the transmission amplitude remains unchanged. It does not depend 
upon the position of the barrier, but only upon its width (L) and height (Vq), apart, of course, upon 
the value of the energy E. 



Now consider the twin barrier problem, the shape of the potential is drawn in Fig. 1-b. The 
particle approach stems from considering the hypothetical case of wave packets such that their spatial 
dimensions are small compared to d. 

Consider an incoming wave packet from the left. When it reaches the first barrier at a; = 0, it will 
produce a first reflected wave in the region a; < and a first transmitted wave into the the inter-barrier 
separation. For each momentum component fc, these are given by the 

R and T 



of the single barrier given above, with consequent conservation of probability. Now the wave packet 
in the inter-barrier region will travel to the second barrier and then be partially reflected by and 
partially transmitted through the second barrier. These amplitudes will be, respectively, 



2ik(L + d) 

TxRe and TxT 



The latter amplitude is the first contribution to the total transmission amplitude. Again probability 
is clearly conserved at this step. The reflected wave from the second barrier now travelling to the left 
impinges upon the first barrier at x = L and, after its transmission through the first barrier, creates 
the second contribution to the total reflection amplitude. 



2ik(L + d) 

TRe xT 



The reflected amplitude (at x = L) which moves back from the flrst to second barrier is 



2ik!L + d) „ -2ikL 

TRe X Re 



Again probability is conserved. Consequently, the second contribution to the total transmission 
amplitude is then given by 



2ikd ^ 

TR e xT 



The procedure repeats continuously and the sum of these contributions gives 

oo ,1 

2 2ik(L + d) X — ^ / 2 2ikd\ 2 2ik(L + d) / / 2 2ikd\ , ^ 

R + RT e zJ ( ^ ^ j =R^RT e / il-R e j , (5) 

(2) 2 X — ^ / 2 2ikd\ 2 / / 2 2ikd\ , , 

= T Y^{r. e ) / (l^R e ) , (6) 

n=0 



(2) 



T 
^ s 



where the bracketed upper script (2) on the l.h.s. indicates the number of identical barriers and the 
subscript S refers to an infinite sum. The second equality in each of the above equations is legitimate 
because \R e | < 1- 

For an incoming wave packet which is small in spatial dimensions compared to the inter-barrier 
potential distance, the individual contributions used in the sums above should not be added because 
they are incoherent. Each will correspond to a separate outgoing (transmitted or reflected) wave 
packet, 

oo 

= ii?r+ii?rirrEi^i"=2|i?i7(i+ii?r) , (?) 

oo 

T,!L. = iTr^ii?r = iTi7(i+ii?f) . (8) 

n=0 

This is what we call the particle limit. Since probability is conserved at each step of the above 
calculation, probability is conserved overall, 

7^''' +r''' =1 

'^particle ^ 'particle ^ 5 



but the transmission and reflection probabilities in the particle limit arc not those of the wave limit 
when each individual particle contribution is completely coherent [15, 17, 18]. In this case, 



n 



(2) 



2 \Rf [1 + cos(2a)] / [l + + 2 cos(2a) 



tJI = ITI" / [l + \Rf + 2 \Rf cos(2a) 



(9) 
(10) 



with 

and, as expected. 



+ kd 



1 



This result must not be misinterpreted as proof of a single outgoing transmission/reflection amplitude. 
The transition between the wave and the particle limit will be discussed in more detail in the next 
section. 

Now let us return to the transmission amplitude given in Eq.(6), seen in the wave limit. 



(2) 2 / 

Ts =T / D , 



(11) 



with 



D ^l-R e 



l + \R\'e 



From this form, one immediately derives the condition for resonance tunnelling for which the barriers 
becomes "transparent" . The maximum of the above expression occurs when the modulus of the 
denominator is a minimum, i.e. when 



cos(2q!) 



-1 cosia) = 



For these values of a 



Tg'' [cosa = 0] = r'/ |rf ^ 7;''i [cos a = 0] = 1 



(12) 



(13) 



Observe that the transmission probability in the particle limit, Eq.(8), can never be unity, since 

(2) 

|T| < 1. We may also determine the minimum value of ^^^^ which occurs when 

cos(2q;) = + 1 ^ cos(a) = ± 1 . 

For these values of a 



Tf' [cosa = ±1] = r'/ (^1 + ^ 7;'l[cosa = ±l] = l/ (^l + 2|i?|V|r|'^' 

III. TRANSITION BETWEEN WAVE AND PARTICLE LIMIT 



(14) 
(15) 



In this section, we intend to study numerically the behavior of a particular incident wave packet 
impinging on a double identical barrier potential (Fig. 1-b) from the left. The results will graphically 
exhibit the transition from the wave to particle limits and also verify the validity of some of the 
expressions derived in the previous section. For the following numerical calculation, we shall express 
all quantities in terms of the barrier height Vq and/or rnVJ,. For our purposes these values need not to 
be explicitly fixed. We shall consider below a set of barrier widths, the incoming wave packet's mean 
momentum (fco)/energy {Eg) and a set of its momentum spread. The inter-barrier distance will be 
left as a continuous variable in our calculations. The incident wave packet is obtained by superposing 
the planes waves 



exp 



kx 

n 



with the coefficient 



(27r) 



3/4 



exp 



(fc-fco)^ 



which correspond to a Gaussian function centered at kg = ^/2rnE^/h multiphed by a numerical factor 
which normahzes the wave function [8]. The transmitted wave packet can be written as 



*t™(a;,i) 



i2nr 



(2) 

dk T<5 exp 



-(fc-fc„) 



k X — 

h 



where w = \/2'mVo/h . Consequently, the probability of transmission is given by 



(16) 



da; \-i'tra{x,t)\ 

2L + d 



(17) 



Now we are interested in the total transmission probability. This is formally the limit of the above 
for t — > oo. However, there is an alternative procedure to obtain the same result. We send the lower 
X limit to — oo. This picks up all the future (at time t) probability contributions, i.e. phantom wave 
packets. An example of this phenomena is seen in numerical calculations with, say, the step potential 
for times before total {E < Vg) wave packet reflection. If one plots the reflected wave amplitude, 
it does not appear yet in the free potential region where destructive interference occurs but a wave 
does appear for x values in the forbidden region. Of course, the reflection coefficient does not apply 
there, and this traveling wave packet is only virtual or phantom. However, it eventually emerges as a 
physical wave packet in the free region at the appropriate reflection time. The above extension of the 
lower X integration limit is useful for Gaussian convolution because it yields a Dirac delta function. 
This allows us to immediately perform one of the momentum integrals. It in turn eliminates the time 
dependence as it must, and yields the following expression for the total transmission probability. 



\/2mVo a/h 



V2tt 
\/2mVo a/h 



hk 



V2mVo 
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exp 
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{\/2mVo a/h) fhk — hk„ 
V^mVT 



{V2mVoa/hf 




(18) 



From the previous equation, observing that is a function of E/Vg, \/2mVo d/h and \/2mVo L/h, 

(2) . 

we immediately see that Vr^ is completely determined once fixed the following four adimensional 
quantities. 



Eo y/2mVoa ^/2mVod 
Vo ' h ' h 



and 



V2mVoL 



is calculated numerically and plotted in Fig. 2, where we have chosen the ratio Eg/Vo equals to 
1/2. In the upper part of the figure. Fig. 2-a, we set the barrier width, 

^/2mVoL/h^ 1 , 



and display the total transmission probability versus the inter-barrier distance, y'2mVo d/h, for various 



T 



is 



wave packet widths. The wave limit is obtained when each individual contribution to 

completely coherent. This happens when the inter-barrier potential distance (d) is small if compared 
to the wave packet spatial dimension (a). 



(2) 



where the subscript indicates that 7^^^^ at E 
resonance structure with maxima at 



Eq. For plane waves, we would have a periodic 
v/2mK d/h = (2n + 1) TT / \/2 . 



As shown in Fig. 2-a, the only close approximation to a resonance/transparency occurs for the first 
maximum and for the largest spreads of wave packet plotted. As ^/2mVo d/h increases in our plot. 



for each choice of wave packet width the probabihty tends to a constant value. This constant value 
is equal to 0.46 and it is in agreement with the particle limit given in Eq.(8). This is due to fact that 
for an incoming wave packet which is small in spatial dimensions (a) compared to the inter-barrier 



potential distance (d) the contributions to 



[a<d] 



are incoherent, 



In Fig. 2-b, we set the wave packet width, 



x/2mVo a/h = 30 , 



and display the total transmission probability versus the inter-barrier distance, y^2mVo d/h, for various 
barrier widths. Again the gradual transition from wave to particle limit is exhibited by the damping 
of the oscillations. Each particle limit is different because the values of R and T depend upon the 
barrier width. Also it is to be noted that the maxima and minima of each curve occur at the same 
values of \/2mVo d/h. This is not true in general. It happens here because of our choice oi Eg = Vo/2. 
This particular ratio results in 0=0. In general </> is dependent upon the barrier width and this would 
separate somewhat the maxima/minima. 

The curves, plotted in Fig. 2, clearly show the transition from the wave to the particle limit. 
However, another even more direct way to show the particle limit is to display the probability density 
as a function of x for the transmitted wave. In Fig. 3, this is shown for a case in which multiple wave 
packets exist and for a time at which two have emerged. These are numerical calculations and we 
take the opportunity to compare the position of these first two maxima with that predicted by the 
stationary phase method, SPM, approximation [19] applied to the particle expression of Eq.(6), 



T R 



This method gives us the times of exit of the maxima and, knowing the group velocity hk„/m, we can 
calculate their later positions. The phase of the first transmitted wave packet is extracted from 



The SPM then gives the position of the first maximum at time t, 



2L 



2^ 
dk 



with the derivative calculated at the maximum of the Gaussian distribution, kg. In this step we have 
neglected the shift in this maximum produced by the transmission amplitude. The derivative term 
is proportional to twice the transition time through a single barrier. It is thus proportional to the 
time needed to tunnel through the twin barriers. Note that this contribution exists even if, as in 
our case, = 0. The SPM transition time for a single barrier is surprising because for very large L 
it becomes independent of the barrier width. This is the Hartman effect [10], which has become of 
renewed interest in recent years [11,12]. 

The phase of the second transmitted wave packet is obtained from 



X R 



whence 



dk 



2d 



The derivative terms here are related to the refiection delay times and, as is well known, are seen to 
be equal to the tunneling times. These SPM position of the maxima, for the case plotted in Fig. 3, 
are indicated by the vertical lines. Agreement with the numerical calculation is excellent. 

However, notice that if the SPM is applied to the summed (wave limit) expression it will yield a 
single maximum which does not coincide with any of the physical particle maxima. The use of the 
SPM must be augmented with, at the very least, a knowledge of the number of maxima involved. 



IV. THREE BARRIER ANALYSIS 



The three identical barrier problem, Fig. l-(c), is somewhat more complicated than the twin barrier 
case, particularly in the extraction of the conditions for resonance/transparency. We have now two 
identical inter-barrier regions and hence different processes or paths can contribute to coherent out- 
going "particles", i.e. with the same exit times. These terms must be summed before squaring the 
amplitudes. For brevity, we derive only the transmission amplitudes. 

We proceed as follows. The leading transmission term is that without any internal reflection, i.e. 
T . After including all internal reflections, we find 



3 r 2 2ikd 

T \ l + 2R e 



4 4ikd 

3R e 



(19) 



_3 _2 2ikd 



The factor two before the T R e term allows for a double reflection, R , in each of the inter-barrier 
regions. The coefficient (three) in the third term allows for a four-fold reflection, R , in each of the 
inter-barrier regions plus a double reflection, R' x i?', in both regions. Using R ^ —i \R\ e"^, Eq.(19) 
can be rewritten as follows 

(20) 



3 / / , ,2 2ia \ ^ 3 / : 

T / (l + \R\ e ) /D 



However, this is by no means all, we can also have reflection from the third barrier, backward tunnelling 
through the second barrier and reflection from the first barrier and again tunnelling through the second 
barrier before exiting the structure. This contribution (to leading order) yields T R e * ^ ' . When 
all additional internal reflections arc allowed for, we obtain the partial sum of 



5 2 2iA;(Z, + 2 

T R e 



d) ( 2 2ikd 

|l + 4i? e 



4 4.ikd 

lOi? e 



which leads to 



5 2 2ifc(L-f2d) 

T R e 



/ 4 3 , ,2 , ,2 iia I 4 

/D = -T \R\ \T\ e /D 



(21) 



(22) 



term, 



The higher denominator power (-D^) follows from the fact that each inter-barrier region is crossed 
twice m the leading T R e path. This process continues for the T R e 

7 4 4ifc(L + 2d) / 6 3 , ,4, ,4 8ia / G 

T R e /D =T \R\ \T\ e /D , 

and so forth. The final result summing Eq.(20), Eq.(22) and Eq.(23) is 
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T 



\R\ \T\ 



\R\ \T\ 
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D 1 



\R\ \T\ 



(24) 



To determinate the resonances, in the wave limit, we first take the modulus squared of the denominator 



D 



\R\ \T\ e 



= 



= H-5|i?| +4|i?| (^1 + |i?| ) cos(2a) -I- 2 cos(4a) 
and then differentiate with respect to a, equating to zero, 

sin(2a) 1 + + 2 cos(2a) 

This results in 

• minimum values of the transmission probability when sin(2a) = 0, which yields 

rf[cosa = 0] = t/\T\\ 

rjl[cosa = o] = |rf 

Ts"[cosa = ±l] = t'/ (^l + 3|i?f) , 
7;|'l[cosa = ±1] = \T\' 



V (l + 4|i?|V|T| 



(25) 
(26) 



• maximum values of the transmission probability when cos(2a) = — ^1 + j /2, which yields the 



resonances 



rli [cosa = ±|T|/2] 



V. CONTINUITY EQUATIONS AND MATRIX METHOD 



(27) 



The procedure for calculating the particle terms for the case of more than three identical barriers 
becomes very difficult. Not only must we be careful of coherence but also calculate enough terms 
to extrapolate the total series. So, in this section we provide an alternative procedure based on the 
matrix method [1,8]. We derive the expression for the wave transmission and reflection amplitudes 
for identical barriers, Fig. l-(d), and then, based upon our results in the previous sections, we 
devise a procedure for deriving the particle sums. We shall then apply this method to the — 4 
transmission amplitude. 

We first define two 2x2 matrices: W[S, x] and the diagonal A [5a;], 



W[5,: 



-Sx 



W[S, 0] 











= W[6, 0] A[(5a 



5e*^ -Se-^"" 

Let the plane wave solutions in the free region to the left of the s-barrier be 
and that within the barrier 

Thus, the plane waves in the extreme left region will be given in terms of Aa (incoming) and A'^ 
(reflected). The final factors on the extreme right of an A^-barricr system will be A„ (transmitted) 
and A'^ (to be set to zero below). 

The continuity equations at the two edges of the s-barrier read 

W[ik,{s-l){L + d)] = W[p,{s-l){L + d)] [B, B'^]' , 

W[p, sL + [s - l)d] [Bs B^]* ^ W[ik, sL + {s- l)d] [A, A'J* , s = 1, 2, ...N . 

Applying successively these matrix equations, we can express the matrix equation between (Aq , A'g) 
and (Ajv , A'„), 



A' 



Y[ W ' [ik, [s -1){L + d)] W[p, (s -1){L + d)] W ' [p, sL + [s - l)d] W[ik, sL + {s - l)d' 

S — 1 

= A[ikd] I A[-ifcd] W'' [ik, 0] W[p, 0] A[-pL] W'' [p, 0] W[ik, 0]|" A[ikNL + ik{N - l)d' 

Introducing the matrix M , 

M = A[-~ikd] W'' [ik, 0] W[p, 0] A[~pL] W'' [p, 0] W[ik, 0] , 
the previous matrix equation can be rewritten as follows 

N 

A[ikd]M A[ikNL + ik{N -l)d] 
A straightforward calculation shows that 
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An 

A'„ 



(28) 



with 



F : 



Te 



.k(L + d) 



and G — 



\T\e' 



T e ' 



\R\ .fcd 



Observe that M has unit determinant, |F| — \G\ ~ 1. For an incoming wave from the left, A'^ = 0, 
Eq.(28) gives 



, ikN(L + d) . 

A„ = ( Af ) e ' ^„ , 



a: 



ik[N L + (N - •2)d ] 

e A. 



(29) 



where the sub-indices specify the matrix element. Thus, the reflection and transmission amplitudes 
are given by 



(N) - 2ikd / N\ I / iV\ 

= e (m ) J{m )^^, 

For two and three barriers, we have respectively 

(m') = + |Gf and (a/) = F (f' + |Gf) + |Gf (F + F*) 
Consequently, in agreement with our previous results. 
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T / D 



and 



1/ \F e 



1 + 



|G| \ , |G| 



F^ 



F* 



These results are easy to find since we already knew the particle expressions. We now derive the 
= 4 transmission amplitude. By using 



M 



F + \G\ ] +|G| (F + F*) , 



from Eq.(30), we get 



Tf . 1/ Fe 



iik{L+d) 



F' 



F' 



(31) 



The leading term will be T / D . This term must be factorized. The remaining expression can and 
must be expressed in terms of \R\ , |r| and e . By adding and subtracting, in the square bracket 
of the denominator of Eq.(31) a D term, and by recalling that \F\ = 1 + |G| , we can rewrite the 
transmission amplitude as follows 



(4) 



/ 4 4ik{L + d 

l/\Fe 
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F \ F 



4 i^k^L + d) 3 \G\' \G\ 
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(32) 



To obtain the particle sum, we must explicit the denominator as a series in the numerator. For A'^ = 4, 
we find up order |r| . 



^ ll-2\R\\"" 



,4 8sQ , ,2 6i, 

R e - R e 



(33) 



For the individual particle contributions we must expand the D terms in the numerator. The above 
result agrees with a direct, but more tedious, particle calculation to this order. The maximum and 
minimum values of the transmission probability are obtained by following the procedure of the previ- 
ous section. This results in 

• minimum values of the transmission probability 



|r|7 (l + 5|i?r/27) , 

i/ (i+8ii?ivirr ~ 



7;J, cosa = ±|T|/V6 
7^^^^ [cos a = ± 1] 
• maximum values of the transmission probability 

7^1, [cos a = 0] = T;*'! [cos a = ± \T\/V2] = 1 



(34) 
(35) 

(36) 



VI. CONCLUSIONS 



We have considered in this paper tunnelling through two, three and four identical barriers. Our 
approach is that of considering the reflection and transmission amplitudes successively for each barrier, 
the so called particle approach. This is directly relevant to situations in which the incoming (single) 
wave packet is small compared to the inter-barrier distance d. However, it also yields the wave limit 
if all particle contributions are summed and considered as a single outgoing reflected and transmitted 
wave. In this wave limit one encounters resonance phenomena. We have re-derived the condition for 
resonances and expressed them in terms of cos a. Resonances do not occur in the particle limit. We 
have exhibited numerically the transition from the wave to particle limit for the case of twin tunnelling 
by plotting the transmission probability as a function of the inter-barrier distance d. In the particle 
limit oscillations die out. 

The calculation of the individual particle terms becomes cumbersome for four barriers and higher. 
Care must be taken to sum coherent contributions which are produced from different paths within the 
potential structure. For these cases it is simpler to invert our procedure and first calculate the wave 
limit amplitude. This has been done, as an example, for the four barrier case in the previous section 
with the help of a matrix method and the rules given for its representation in "particle" form, i.e. in 
terms of R and T and powers of e .A series expansion of the denominators completes the procedure. 



There are a number of consequence of our analysis which merit a mention: 



1) The resonance conditions for multiple identical barrier tunnelling will automatically extrapolate 
into the case oi E > V„, i.e. for above barrier diffusion. Thus, there are two classes of above barrier 
resonances. 



1-a) When above barrier resonance occurs for a single barrier, |r| = 1, independently of the 



(2) ^ 

value of d we have = T , where T is now the above barrier expression, i.e. p ^ iq with 
q = ^y2m{E -V„)/h in Eq.(2). 



1-b) When |T| 7^ 1 but (j) + kd = {n + 1/2) tt as in the tunnelling case. These are additional 
resonances corresponding to constructive interference effects of the twin barriers possible in the wave 



limit even if each single barrier is not resonant. 



2) For twin barriers, the conditions for tunnel resonance requires identical barriers. For two different 
barriers, either in width or height or both, the formula for the transmission amplitude becomes 

with the indices indicating the individual barrier amplitudes. Now, for no value of kd will this be of 
unit modulus. We expect the same to happen also for higher numbers of barriers. 

(2) 2 

3) Consider twin tunnelling at resonance. The phase of is totally given by the phase of T . The 
time calculated at a; = 2L + d, the exit point for transmission, is 



k—ko 

where Vg is the group velocity in the free space. Thus, due to the linear term in d, there is no genera- 
lized Hartman effect [15] at resonance. 

Finally, we wish to place our results in a broader context. The existence of wave and particle limits is 
related to the relative sizes of the incoming wave packet and the size of the potential structure. The 
results of this paper are an example of the physical relevance of wave packet dimensions. Interest- 
ing applications occurs in particle oscillation phenomena [20-23], relativistic tunneling [24-26], laser 
interaction with dielectric blocks [27,28], and, generally speaking, in all interference based results. 
These features are generally neglected in the literature where only single plane waves are used. Plane 
waves are a legitimate operational tool, but working only with single plane waves (infinite wave packet 
size) may obfuscate significant physical insight. 
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Figure 1: Potential shapes: (a) single barrier, (b) twin barriers, (c) triple identical barrier and (d) 
general structure for matrix calculation. In (d) the A/ A' and B / B' terms indicate the amplitudes in 
the free and potential regions. 
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Figure 2: Transmission probability vs the inter-barrier distance. Upper curves are for a fixed barrier 
width (L) and various wave packet widths (a). Lower curves are for a fixed wave packet width and 
various barrier widths. The curves show the gradual transition from the wave limit (continuous line) 
to the particle limit (dotted line) exhibited by the damping of the oscillations. Observe that the 
particle limit is the same for all the curves plotted in Fig.2-a whereas each particle limit is different 
in Fig.2-b. This is due to the fact that T^^^i^i^ depends on the barrier width. 
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Fig. 3 
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Figure 3: The density probability for the transmitted wave is plotted as a function of x for a fixed 
time. This time is such that two outgoing waves are seen. The SPM estimates (vertical lines) are in 
excellent agreement with the numerical calculation. 



